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We analyze the recent results of the MiniBooNE short-baseline experiment on £> M — > i> e oscillations 
in a minimal model-independent framework of antineutrino mixing in conjunction with the positive 
LSND signal and the negative KARMEN measurements. We show that the data of the three short- 
baseline P M — t v £ experiments are compatible. Taking into account also the model-independent 
constraints due to the limits on short-baseline disappearance obtained in reactor antineutrino 
experiments, we find that the favored region of the effective oscillation parameters lies within 2 x 
10~ 3 < sin 2 20 < 5 x 10~ 2 and 0.2 < Am 2 < 2eV 2 . 
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I. INTRODUCTION 

The MiniBooNE collaboration [l| recently reported the 
observation of a signal of short-baseline Cy, — > v e transi- 
tions compatible with that observed in the LSND ex- 
periment [2[ . The agreement of the MiniBooNE and 
LSND signals in favor of neutrino oscillations is remark- 
able, because the two experiments observed the signal 
of 9^ — > v e transitions at different source-detector dis- 
tances and different neutrino energy ranges. Since only 
the ratio of distance and energy is similar in the two ex- 
periments and neutrino oscillations depend just on this 
ratio (see Refs. [H-Qil), the neutrino oscillation expla- 
nation of the two signals is strongly favored. On the 
other hand, the MiniBooNE collaboration did not ob- 
serve any signal of short-baseline — > v e transitions 
11] compatible with the MiniBooNE and LSND signals 
of Vp — > v e transitions. Therefore, it is possible that 
the effective parameters which govern neutrino and an- 
tineutrino oscillations are differen t, may be because of a 
violation of the CPT symmetry [12h30| |. From a phe- 
nomenological point of view, it is interesting to consider 
the neutrino and antineutrino sectors independently, es- 
pecially in view of possible experimental checks of the 
short- baseline P M — > v e signal |3ll - t34| . In this paper we 
adopt this point of view and we present the results of a 
combined fit of the MiniBooNE and LSND antineutrino 
data in favor of short-baseline — > v e transitions, to- 
gether with the constraints imposed by the data of the 
KARMEN experiment [35| in which the transitions have 
not been observed. We also take into account the model- 
independent constraints imposed by the data of reactor 
v e disappearance experiments. 
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In the analysis of the data of — > v e oscillation exper- 
iments we consider the simplest case of an effective two- 
neutrino-like short-baseline oscillation probability, simi- 
lar to that obtained in the case of four-neutrino mixing 
(see Refs. 0,1,1,0), 

P^ e (L/£) = sm 2 2tfsm 2 (^) , (1) 

where Am 2 is the relevant neutrino squared-mass dif- 
ference and § is the effective mixing angle for — > D e 
transitions. 

The plan of the paper is as follows. In Sections |TT] 
and IIIII we present, respectively, the results of the fits 
of MiniBooNE and LSND antineutrino data, and in Sec- 
tion IIVI we discuss the results of the combined fit. In 
Section |V] we present the results of the fit of KARMEN 
data and in Section IVTl we discuss the results of the com- 
bined fit of MiniBooNE, LSND and KARMEN data. In 
Section lVlIl we discuss the implications of the constraints 
from reactor v e disappearance experiments. Finally, in 
Section IVlIII we draw the conclusions. 

II. MINIBOONE 

The MiniBooNE collaboration presented recently the 
results of a search for — > v e oscillations obtained with 
a data sample corresponding to 5.66 x 10 20 protons on 
target [l|. The MiniBooNE detector is located at a dis- 
tance of 541 m from the neutrino source. The neutrino 
energy spectrum for the oscillation analysis ranges from 
475 MeV to 3 GeV. Hence, the ratio L/E from which 
the oscillation probability in Eq. (|T|) depends ranges from 
0.18 to 1.14 m/MeV, leading to a sensitivity to — s- v e 
transitions for Am 2 > 1CP 1 eV 2 appropriate for check- 
ing the signal observed in the LSND experiment 0] (see 
Section Hn|. 

The excess of f> e -like events found by the MiniBooNE 
collaboration agrees with the excess found in the LSND 



2 







MB 


LS 


MB+LS 


KA 


MB+LS+KA 


Re 


(MB+LS+KA)+Re 


No Osc. 




21.4 


15.0 




6.8 




51.0 






NDF 


16 


5 




8 




56 






GoF 


0.16 


0.010 




0.55 




0.66 




Osc. 


Xmin 


11.7 


1.4 


14.6 


6.4 


25.7 


48.5 


77.3 




NDF 


14 


2 


18 


6 


26 


54 


82 




GoF 


0.63 


0.51 


0.69 


0.38 


0.48 


0.69 


0.63 




sm 2v bt 


0.91 


0.0058 


0.006 


0.0010 


1.00 


0.042 


0.014 




Am 2 , 


0.071 


8.13 
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1.50 

2 
0.47 




6.32 

4 
0.18 




3.02 

2 
0.22 



TABLE I. Values of x 2 5 number of degrees of freedom (NDF), goodness-of-fit (GoF) and best-fit values sin 2 2i9bf, Am^ of the 
oscillation parameters obtained from the fit of various combinations of MiniBooNE (MB), LSND (LS), KARMEN (KA) and 
reactor Bugey and Chooz (Re) antineutrino data. The first three lines correspond to the case of no oscillations (No Osc). The 
following five lines correspond to the case — s> oscillations (Osc). The last three lines give the parameter goodness-of-fit 
(PG) [3g]. The variations of sin 2 2i?bf and Am 2 f depending on the fitted data sets are due to the oscillating character of P^ tl ^ i 7 c 
in Eq. ©. 
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FIG. 1. Allowed regions in the sin 2 2§-Am 2 plane and 
marginal A% 2 = X 2 — Xmin' s f° r sin 2 2i? and Am 2 obtained 
from the fit of MiniBooNE (MB) antineutrino data in the en- 
ergy range E > 475 MeV [l[. The best-fit point is indicated 
by a cross. 



experiment at different source-detector distance and neu- 
trino energy, but similar ratio L/E. This is a strong 
indication in favor of the neutrino oscillation explana- 
tion of the two signals, analogous to the confirmation 
of solar neutrino oscillations by the very-long-baseline 
KamLAND reactor experiment [3 71 ] and the confirmation 
of atmospheric neutrino oscillations by the long-baseline 
K2K 38] and MINOS [H| accelerator experiments. 



In this paper we fit the MiniBooNE data reported in 
Fig. 1 of Rcf. lQ using the method and data given in 
the MiniBooNE data release in Ref . [4(| , which are rela- 
tive to the previous MiniBooNE publication [4l[ on the 
search for — > D e (the MiniBooNE data release relative 
to Ref. [U is still not available). We rescaled the signal 
predicted with the method described in Ref. [40] from 
the 3.39 x 10 20 _protons on target corresponding to the 
sample in Ref. [41] to the 5.66 x 10 20 protons on target 
corresponding to the sample in Ref. jlj. The fractional 
covariance matrix of systematic uncertainties should be 
similar in the two data releases. We corrected the statis- 
tical part of the covariance matrix by taking into account 
the different number of background events. In the fit we 
consider not only the P e MiniBooNE data, but also the 
data, which are important because of the correlated 
uncertainties. Since the data obtained with 5.66 x 10 20 
protons on target are not available, we consider th e v a 
data given in the MiniBooNE data release in Ref. [40] ]. 
The difference is not crucial, because the P M data have 
only an indirect effect on the measurement of the — > v e 
signal through the correlated uncertainties. 

The results of the least-squares fit of MiniBooNE data 
are presented in the first column of Tab. U and in Fig. [T] 
The best-fit values of the oscillation parameters and the 
allowed regions in the sin 2 2??-Am 2 plane are similar to 
those obtained by the MiniBooNE collaboration |l| . The 
goodness-of-fit in the case of no oscillations may seem too 
high in comparison with that given in Ref. [lj and not 
sufficient to require oscillations. Since in our calculation 
we fit both the v e and data we have 16 degrees of 
freedom, with \ 2 = 21.4. However, most of the \ 2 is due 
to the i> e data, which have only 8 degrees of freedom, 



1 We would like to thank W.C. Louis for sending us the precise 
values of the data in Fig. 1 of Ref. 
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corresponding to the 8 energy bins in Fig. 1 of Ref. |l| 
for E > 475 MeV. If we restrict the x 2 to the six energy 
bins from 475 to 1300 MeV we have x 2 = 16-8, which 
is similar to the x 2 = 18.5 reported in Ref. [lfl for the 
energy range from 475 to 1250 MeV. Therefore, we agree 
with Ref. |l| on the opinion that a background-only fit 
of MiniBooNE data is disfavored. 

The first column of Tab.|T]shows that the oscillation hy- 
pothesis fits the MiniBooNE data with a x 2 much lower 
than in the case of no oscillations, improving significantly 
the goodness-of-fit. The decrease of the x 2 with respect 
to the case of no oscillations is mainly due to the im- 
proved fit of the six T> e energy bins from 475 to 1300 MeV 
which give a contribution to Xmin 01 7-2,in approximate 
agreement with the 8.0 reported in Ref. [l| for the energy 
range from 475 to 1250 MeV. 

From Fig. [T] one can see that, although the best-fit 
value of sin 2 2d is close to unity, in practice all the al- 
lowed straight region in the log-log plot ranging from 
sin 2 20 « 1 and Am 2 « 6 x 10~ 2 eV 2 to sin 2 2d « 2 x 10~ 3 
and Am 2 « 1 eV 2 , as well as a small area around 
sin 2 2d ~ 7 x 10 -3 and Am 2 5eV 2 , are equally plausi- 
ble. This is important, because large values of the effec- 
tive mixing angle are excluded by the limits on P e disap- 
pearance obtained in reactor antincutrino experiments, 
as explained in Section IVIII 
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FIG. 2. Allowed regions in the sin 2 2i9-Am 2 plane and 
marginal Ax 2 's for sin 2 2$ and Am 2 obtained from the fit 
of LSND antineutrino data. The best-fit point is indicated by 
a cross. 



III. LSND 

The LSND experiment [2J observed an excess of P e 
events coming from possible P^ — > v e transitions in a 
beam of P M produced by /i + decay at rest, 

fJ- + -> e+ + v e + P M . (2) 

The energy spectrum of P M is given by (see Ref. [42]) 

^ tt (E)ocE 2 (3-4E/m^ , (3) 

for E smaller than 

£max - (m„ - m e )/2 ~ 52.6MeV . (4) 



The v e 's produced by Vu 



v e transitions were de- 



tected at an average distance of 30 m through the inverse 
neutron decay process 



v e + p — » n + 



(5) 



in a detector consisting of an approximately cylindrical 
tank 8.3 m long by 5.7 m in diameter, filled with liq- 
uid scintillator. Hence, the oscillation distance varied 
between L m - m — 25.85 m and L max = 34.15 m. The 
cross section of the detection process in Eq. flS} is (see 
Refs. 0,111,111) 



o-p eP (E e ) oc E e p e 



(6) 



where E e and p e are, respectively, the positron energy 
and momentum. Neglecting the small recoil energy of the 



neutron, the positron energy E e is related to the neutrino 
energy E by 



E,, = E 



m n = E - 1.293 MeV, 



(7) 



where m p and m n are, respectively, the proton and neu- 
tron masses. The neutrino energy threshold is given by 



E„ 



(m r , 



2 m z 



1.806 MeV. 



(8) 



The LSND detector had a positron energy resolution 
which we assume to be Gaussian: 



R(E e i Ee 

with [3 



exp 



2*1 



6 Ee = 3.3 MeV 



E„ 



50 MeV 



(9) 



(10) 



We fit the LSND data in Fig. 16 of Ref. 0, which 
gives the measured P e events in five bins of mea- 
sured positron energy, their uncertainties SN^ and the 
expected number of background events B l 9 . The fit is 
obtained by minimizing the least-square function 



5 

E 
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(11) 
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The number N„ of v,, — > v„ events in the ith bin 



Ni _ 



where 

= 1.29 x 10 4 (13) 

is the number of events expected for 100% — > v e 
transmutation in the measured energy range between 
Ef n = 20 MeV and £™ ax = 60 MeV. We obtained this 
number by dividing the — » v e excess in Fig. 16 of 
Ref. y| (34.1 events) by the probability in Tab. XI of 
Ref. @ (2.64 x 10~ 3 ). The factor rj in Eq. (JTTJ is intro- 
duced in order to take into account the relative uncer- 
tainty of 0: 

5r] = 0.1 . (14) 

Figure [2] shows the allowed regions in the sin 2 2i9- Am 2 
plane that we obtained from the minimization of % 2 , 
with the best-fit values of the oscillation parameters in 
the second column of Tab. [J] The allowed regions in 
Fig. [5] are similar to those presented by the LSND col- 
laboration in Ref. [2|, with some differences due to the 
fact that we fitted a data set which is smaller than that 
used by the LSND collaboration. Similar problems have 
been encountered in the fits presented in Refs. 0, |45j . 
In particular, the allowed straight region in the log-log 
plot ranging from sin 2 2?9 rj 1 and Am 2 rj 5 x 10 _2 eV 2 
to sin 2 2d rj 10~ 3 and Am 2 rj 2 eV 2 and the allowed 
region at large values of Am 2 are similar to that ob- 
tained by the LSND collaboration (see Fig. 27 of Ref. 0). 
There is some discrepancy at intermediate values of Am 2 , 
where we find two favorite regions at Am 2 rj 5 eV 2 and 
Am 2 sa 8eV 2 , which however are similar to those ob- 
tained in Ref. [45] from a fit similar to ours. Hence, we 
think that the allowed regions in Fig. [2] are fairly repre- 
sentative of the parameter space allowed by the LSND 
data. 



IV. COMBINED FIT OF MINIBOONE AND 
LSND DATA 

Comparing Figs. Q] and [2] one can see that there is a 
remarkable agreement between the MiniBooNE and the 
LSND allowed regions in the sin 2 2i?-Am 2 plane. The 
results of the combined fit are given in Fig. [3] and in 
the third column of Tab. |U The excellent parameter 
goodness-of-fit (PG) [36[ of the combined fit quantifies 
the good compatibility of MiniBooNE and LSND data. 



of measured positron energy between E\ — AE e /2 and 
El + AE e /2 is given by 
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FIG. 3. Allowed regions in the sin 2 2i9-Am 2 plane and 
marginal Ax 2 's for sin 2 2i3 and Am 2 obtained from the com- 
bined fit of MiniBooNE (MB) and LSND antineutrino data. 
The best-fit point is indicated by a cross. 



From Fig. [3] one can see that the combined fit favors the 
allowed straight region in the log-log plot ranging from 
sin 2 2i? rj 1 and Am 2 rj 5 x 10~ 2 eV 2 to sin 2 2d rj 2 x 10~ 3 
and Am 2 rj 1 eV 2 and an island at sin 2 2d rj 6 x 10~ 3 
and Am 2 rj 5 eV 2 . 



V. KARMEN 

The KARMEN experiment [HJ searched for — > v e 
transitions using a beam of produced by the process 
of /i + decay at rest in Eq. @ and detected through the 
inverse neutron decay process in Eq. ([5]). Since these 
processes are the same as those in the LSND experi- 
ments, the fit of KARMEN is analogous of that described 
in Section [m] for the fit of LSND data. In the KAR- 
MEN experiment the oscillation distance varied between 
^min = 15.935 m and L max = 19.465 m. The energy res- 
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olution uncertainty was 



5 Ec = 0.115 MeV 



MeV 



(15) 



We fit the KARMEN data in Fig. lib of Ref. [35J, which 
gives the measured v e events in nine bins of measured 
positron energy and the expected number of background 
events . The number of events expected for 100% 
D„ — > v e transmutation was 



5826(1 ±0.0923) 



(16) 



FIG. 4. Exclusion curves in the sin 2 2$-Am 2 plane obtained 
from the fit of KARMEN antineutrino data. 



Since in some energy bins in Fig. lib of Ref. [35j the 
number of measured events is zero, we perform the fit by 
minimizing the least-square function [46J 



= 2E 



Bl + r,Nl -m +m In 



(17) 



with 5r) = 0.0923 from Eq. (TT6l) . 

Figure [4] and the fourth column of Tab. U give the result 
of the fit of KARMEN data. The best-fit values of the 
oscillation parameters in Tab. U and the exclusion curves 
in Fig. @] are in agreement with those found by the KAR- 
MEN collaboration [35II. a s well as with the results of the 
fits presented in Refs. [44l. |45| . 



VI. COMBINED FIT OF MINIBOONE, LSND 
AND KARMEN DATA 

The results of the combined Fit of MiniBooNE, LSND 
and KARMEN data on Dn — > v e oscillations are given in 
the fifth column of Tab. U and in Fig. [5] 

Comparing Figs. [3] and [5] one can see that the inclu- 
sion in the analysis of KARMEN data has mainly the 
effect of disfavoring the regions at Am 2 > 3 eV 2 allowed 
by MiniBooNE and LSND data. The straight region in 
the log-log plot ranging from sin 2 2$ « 1 and Am 2 rj 
5 x 10~ 2 eV 2 to sin 2 2tf « 10~ 3 and Am 2 rj 2 eV 2 allowed 
by MiniBooNE and LSND data suffers only a small push 
towards smaller values of sin 2 2*&. The best-fit point lies 
in this region, close to the large-sin 2 2$ and small-Am 2 
edge. However, from the marginal A\ 2 = X 2 — Xmin' s 



for sin 2 2d and Am 2 one can see that in practice all the 
straight region is equally favored. This is important for 
the compatibility with the reactor limits on sin 2 2i9 dis- 
cussed in the next Section. 



VII. CONSTRAINTS FROM REACTOR v e 
DISAPPEARANCE 

Several reactor experiments have searched for the 
short-baseline disappearance of z/ e 's (see Refs. [Tol. l43l|). 
without positive results (apart from the hint discussed in 
Ref. 03)- Such a lack of short-baseline v e disappearance 
constrains the probability of all transitions of v e to other 
flavor antineutrinos and all transitions from other flavor 
antineutrinos to v e . We are interested in particular in 
— > v e oscillations, which are of the second type. The 
constraint is model-independent and does not require any 
assumption on the type of mixing and on the number of 
massive neutrinos, because it follows from simple particle 
conservation, which is a characteristic of oscillations. In 
fact, since in neutrino oscillations a v e must come from 
an antineutrino of some flavor, the sum over the proba- 
bilities of transition of any flavor antineutrino into v e is 
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FIG. 5. Allowed regions in the sin 2 2i?-Am 2 plane and 
marginal Ax 2 's for sin 2 2i9 and Am 2 obtained from the com- 
bined fit of MiniBooNE (MB), LSND and KARMEN antineu- 
trino data. The best-fit point is indicated by a cross. 



equal to unity: 



* . = i 



Then we have the inequality 



(18) 



(19) 



Hence the lower limits obtained in short-baseline reac- 
tor antineutrino experiments on Pp e _>p e imply model- 
independent upper limits on P B _>.p e . 

Considering the simplest case of an effective two- 
neutrino-like short-baseline v e survival probability which 
is governed by the same Am 2 relevant for the effec- 
tive short-baseline probability of — > v e transitions in 
Eq. ([T]), we have 

P^ c {LjE) = 1 - sin 2 2d ee sin 2 {^^\ , (20) 

where d ee is the effective mixing angle, which can be 
different from that of £>„ — > v e transitions (which we 
have denoted for simplicity d, but could have been called 
more appropriately d efJl ). In this case, the inequality in 
Eq. CEH]) implies 



sin 2 2d < sin 2 2d, 



(21) 



Therefore, the exclusion curves obtained in short-baseline 
reactor antineutrino experiments which place upper lim- 
its on the value sin 2 2d ee as a function of Am 2 imply 



FIG. 6. Superposition of the allowed regions in the sin 2 2i9- 
Am 2 plane obtained from the combined fit of MiniBooNE 
(MB), LSND (LS) and KARMEN (KA) v„ -> v e data and 
the exclusion curves obtained from the fit of reactor Bugey 
(Bu) and Chooz (Ch) v a — >• u £ data. 



model-independent upper limits on the value of sin 2$ 
in short-baseline — > v e experiments. 

Figure [5] shows a superposition of the 90% and 99% 
C.L. allowed regions in the sin 2 2d- Am 2 plane obtained 
from the combined fit of MiniBooNE, LSND and KAR- 
D e data and the exclusion curves obtained 



MEN Dn- 

in Ref. 47J from the fit of reactor Bugey [48[ and Chooz 
49] v e -¥ P e data, which currently provide the most strin- 
gent constraints on short-baseline reactor D e disappear- 
ance. The inequality (f2"Tj) implies that in Fig. [5] the large- 
sin 2 2d part of the straight region below Am 2 s» 2 eV 2 
allowed by the combined fit of MiniBooNE, LSND and 
KARMEN — > v e data is excluded by the results of re- 
actor antineutrino experiments. Quantitatively, only the 
parts with sin 2 2i9 < 3 x 1(T 2 and sin 2 2d < 5 x 10~ 2 are 
allowed at 90% and 99% C.L., respectively. 

The inequality (|2"Tj) constrains the effective amplitude 



sin" 2d of short-baseline P M — > D e transitions, but does 
not allow a combined fit of accelerator — ¥ v e data and 
reactor v e — ¥ D e data. Such a combined fit can be done 
if the upper limit in Eq. (|2ip applies, i.e. if the inequal- 
ity (j2"Tj) effectively becomes an equality. This is the case 
if Pv^ve 3> -Pp q ^p c f°r oc ^ e, fi. In the following we 
consider this interesting possibility, which allows us to 
combine the accelerator and reactor data in order to find 
the preferred region in the space of the oscillation pa- 
rameters which could be explored by future experiments 

sis 
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FIG. 7. Allowed regions in the sin 2 2i9-Am 2 plane and 
marginal Ax 2 's for sin 2 2$ and Am 2 obtained from the com- 
bined fit of MiniBooNE (MB), LSND (LS) and KARMEN 
(KA) — y v e data and the exclusion curves obtained from 
the fit of reactor Bugey and Chooz (Re) v e — > v e data. The 
best-fit point is indicated by a cross. 

Figure [7] and the last column of Tab. U give the results 
of the combined fit of accelerator MiniBooNE, LSND and 
KARMEN — > D e data and reactor Bugey and Chooz 
v e — ¥ v e data assuming an equality in Eq. (|21l) . The value 
of the parameter goodness-of-fit in Tab. U shows that the 
accelerator and reactor data are compatible under the 
hypothesis of —¥ v e oscillations. 



From Figure[7]one can see that there is a favorite region 
at about 95% C.L. around the best-fit point for 2x 10~ 3 < 
sin 2 2d < 5 x 1CT 2 and 0.2 < Am 2 < 2 eV 2 . Larger values 
of Am 2 are allowed only at more than about 95% C.L. 
for 7 x 1CT 4 < sin 2 2d < 5 X 1CT 3 . 



VIII. CONCLUSIONS 

We have considered the recent results of the Mini- 
BooNE experiment [l| on short-baseline — > D e oscilla- 
tions, which confirm the positive LSND signal @- Con- 
sidering the simplest case of an effective two-neutrino- 
like short-baseline oscillation probability which depends 
on only two effective oscillation parameters, sin 2 2$ and 
Am 2 , we performed a combined fit of MiniBooNE and 
LSND data in order to find the allowed regions in the 
parameter space. 

We considered also the results of the KARMEN exper- 
iment 35], in which — > D e transitions have not been 
observed. We have shown that the combined fit of Mini- 
BooNE, LSND and KARMEN data is acceptable and 
leads to a shift of the region allowed by MiniBooNE and 
LSND towards small values of sin 2 2d. 

Finally, we have considered the model-independent 
bound on short-baseline — > v e implied by the lim- 
its on short-baseline v e disappearance obtained in re- 
actor experiments. From a combined fit of accelera- 
tor — > v e data and reactor v e — > v e data we have 
found that, if the — > D e channel is dominant over 
other channels of flavor transitions into v ei the favored 
region of the effective oscillation parameters lies within 
2 x < sin 2 2d < 5 x 1CT 2 and 0.2 < Am 2 < 2eV 2 . 

This region is interesting for a study of the possibilities 
to check the LSND and MiniBooNE indication of short- 
baseline — > v e oscillations with future experiments 
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